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ABSTRACT 

An information-geometric approach to sensor management is in- 
troduced that is based on following geodesic curves in a manifold 
of possible sensor configurations. This perspective arises by ob- 
serving that, given a parameter estimation problem to be addressed 
through management of sensor assets, any particular sensor config- 
uration corresponds to a Riemannian metric on the parameter man- 
ifold. With this perspective, managing sensors involves navigation 
on the space of all Riemannian metrics on the parameter manifold, 
which is itself a Riemannian manifold. Existing work assumes the 
metric on the parameter manifold is one that, in statistical terms, cor- 
responds to a Jeffreys prior on the parameter to be estimated. It is 
observed that informative priors, as arise in sensor management, can 
also be accommodated. Given an initial sensor configuration, the tra- 
jectory along which to move in sensor configuration space to gather 
most information is seen to be locally defined by the geodesic struc- 
ture of this manifold. Further, divergences based on Fisher and Shan- 
non information lead to the same Riemannian metric and geodesies. 

Index Terms — Information geometry; Sensor management 

1. INTRODUCTION 

The work of Amari and others (T) on the use of methods of Rie- 
mannian geometry to analyze statistical estimation problems is of 
increasing interest to researchers in signal processing. This method- 
ology, known as information geometry, provides a rigorous frame- 
work for measuring the power of data to discriminate values of pa- 
rameters. These ideas date back to Rao |9|, who showed that the 
Fisher information of a likelihood used in an estimation problem can 
be seen as a Riemannian metric on the parameter manifold. 

This paper brings an information-geometric perspective to a class 
of sensor management problems by casting the objective of sensor 
management as parameter estimation and describing how this leads 
to the role of sensor management as selecting a Riemannian metric 
for the parameter manifold. Established results in Riemannian ge- 
ometry |3|, outside the context of information geometry, show that 
the collection of all Riemannian metrics on a Riemannian manifold 
is itself an (infinite-dimensional) Riemannian manifold. In problems 
where the collection of possible sensor actions is suitably modeled 
by a smooth finite-dimensional manifold, the space of interest is 
a finite-dimensional sub-manifold of this infinite-dimensional Rie- 
mannian manifold. A perspective is developed in which the best 
sensor management action to take, in terms of gathering the most 
information relevant to the estimation objective, is locally charac- 
terised in terms of geodesic curves in this space. 

Much of the development in subsequent sections of this paper is 
rather abstract and draws upon mathematical machinery that is unfa- 
miliar to many researchers in sensor management area. To provide a 



more concrete context in which to illustrate some of the concepts that 
arise in later sections, is it helpful to begin by setting forth an exam- 
ple problem. Suppose two mobile sensor platforms and one station- 
ary target (emitter) are located in the plane K 2 , as depicted in Fig. 1. 
The goal is to estimate the position of the target from bearings-only 
measurements taken at the sensors. Since the sensors are mobile, the 
sensor management problem is to identify the trajectories of sensor 
motion that will yield the best estimate of the target position. More 
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Fig. 1. An illustrative scenario involves estimating the position 
(x e ,y e ) of a single stationary emitter from bearings-only measure- 
ments received at two mobile sensing platforms located at respective 
coordinates (xi,yi) and (x2, 2/2) in the plane. 

specifically, the target position is (x e ,y e ) and the sensor positions 
are (xj, yj) for j = 1, 2. Denoting Xj = Xj — x e and y"j — yj — y e , 
the bearing of the target from sensor j is ipj = arctan(yj/x :) ). The 
sensor measurements are independent and von Mises distributed, 
each with common inverse dispersion parameter n and with the cir- 
cular mean of the measurement at sensor j having circular mean ipj . 

The following sections proceed to describe the nature of a sen- 
sor model from an information-geometric viewpoint, to define the 
parameter manifold, the sensor manifold, its metric structure, and 
to derive a differential equation that characterizes geodesic curves 
on the sensor manifold. The development departs from the purely 
geometric treatment in 1 3 1 in that it allows for informative prior dis- 
tributions on the parameter manifold rather than restricting attention 
to a volume form corresponding to the Jeffreys prior. Further, the 
Riemannian metric with respect to which geodesies maximize "en- 
ergy integrals" on the sensor manifold is shown to arise from both 
Kullback-Leibler and mutual information perspectives. Throughout 
this process, the example just introduced will be used to illustrate 
these concepts in a concrete fashion. 



2. SENSOR MODEL 

Consider the problem of estimating a parameter 8 from data x col- 
lected by sensors up to time t. Beginning with a prior probability 
distribution for 8, which may reflect what is known from previous 
measurements or side information, the effect of taking measurement 
at time t is to provide a posterior probability distribution, which 
will be assumed to be represented by a posterior probability den- 
sity p(8\x). If the option exists to use one of a parametrized set of 
sensors or sensor configurations, each of these will produce its own 
posterior density. When these posteriors are known, selecting a sen- 
sor configuration to use for a measurement amounts to choosing one 
of them from which to estimate 8. The parameter 8 will henceforth 
be assumed to be an element of an m-dimensional smooth manifold 
M (C°° will be assumed, although C 2 is sufficient for most of the 
discussion here), which will be called the "parameter manifold." In 
the example problem, M is R 2 because the parameter 8 = (x e ,y e ) 
is a physical location in the x — y plane. 

Denoting by I — log p(x\9) the log-likelihood for this problem, 
the Fisher information is 

Fe = Ep(.\ B )[del® del], (1) 

where dgl denotes the derivative of £ with respect to the parameter 
8. This is well known to be equivalently expressed as 

Fe = -E pm [V 2 e i\. 

In this expression, V# represents the covariant derivative along any 
connection in M; all choices of connection give the same quantity 
as {T}. Fg is always a non-negative definite m x m matrix, and in 
what follows it will be assumed to be non-singular, thus giving rise 
to a continuous family of inner products over the tangent spaces of 
the manifold. Direct calculation in the example problem shows that, 
in the coordinate system described above and depicted in Fig. 1, Fg 
has the form 

if) 

= kA ( k ) 4l (&< ® (Vh-Zj) ( 2 ) 
= kA(k) -— j (sin (pj, — cos (fij) ® (siiKfj, — cos <pj), 

where i? 2 = x] + y]. 

Through this mechanism, the choice of a particular sensor leads 
to the association of a positive definite matrix with each 8 £ M, 
thereby imbuing M with a Riemannian metric that measures the 
ability of that sensor's data, at least locally, to discriminate between 
parameter values. It is possible to calculate the shortest distance, in 
terms of this metric, between two values 8. As discussed in [ 1|, the 
Kullback-Leibler divergence between p(x\8) and p(x\8 r ) is approx- 
imately half of the square of the distance between 8 and 9' when 8 
and 9' are close. 



in (3). A point in A4 is a Riemannian metric on M; i.e., it asso- 
ciates a positive definite form gg with each 8 £ M. Under suitable 
assumptions, a metric on M is defined by 

G g {h,k)= Tr(gg 1 h e gg 1 kg) vol(g g ), (3) 

J M 

where vol(<7e) = \J det(<?s) d8. Specific assumptions guaranteeing 
finiteness of this integral are beyond the scope of this discussion, and 
it will suffice for the purposes here to assume directly that it is finite. 

Although the nature of M appears formidable, realistic sensor 
management problems do not require one to work with this entire 
space, rather with a finite-dimensional sub-manifold that inherits the 
metric l|3} from A4 . The assumption that leads to this situation is that 
the collection of all possible sensor configurations is parametrized by 
a smooth manifold S, which will be called the "sensor manifold." In 
the example problem, the sensor configuration is completely speci- 
fied by the positions of the two sensor platforms in the plane; i.e., 
by a — (xi,yi,X2,y2) £ R 4 - In this case, the sensor manifold 
is S = R 4 and the only elements of M of relevance are those 
metrics on M that arise from a sensor configuration a in this four- 
dimensional manifold. 

Beginning with a sensor configuration a G S gives rise first to 
a likelihood p a (x\8) and consequently to a Riemannian metric g(a) 
on the parameter manifold M, as described in Section|2] As a Rie- 
mannian metric on M, g(a) is an element of A4(M). This mapping 
g : S — > M taking a to g(a) will be called the "sensor geometry." 
In what follows, g will be assumed to be smooth and one-to-one and 
g(S) a sub-manifold of M. Weaker assumptions are possible (e.g., 
g is an immersion), but full generality is is not needed here to ade- 
quately illustrate the method. Through the sensor geometry map, the 
finite-dimensional manifold S inherits the Riemannian structure of 
M; i.e., the distance between two sensor configurations in a\ and 
o"2 in S is taken to be the distance between g(p\) and g(p%) in M. 
This construction endows the sensor manifold S with its own Rie- 
mannian metric which captures, in information-theoretic terms, the 
"complementariness" of sensor configurations. 
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Fig. 2. The sensor geometry map g allows the finite-dimension man- 
ifold S of sensor configurations to inherit a Riemannian metric from 
the infinite-dimensional manifold A4 of all Riemannian metrics on 
the parameter manifold M. 



3. THE SENSOR MANIFOLD 

It has been shown 1 3 1 that the collection M (M) of all Riemannian 
metrics on the manifold M is an infinite-dimensional (weak) Rie- 
mannian manifold. The structure of its tangent space is described 



4. GEODESICS 

The objective of determining good trajectories in the sensor mani- 
fold S will be addressed by relating these to geodesic curves in M. 
Following (3), consider a smooth curve 7 : [0, 1] — > M. For each 



t G [0, 1], 7(f) is a Riemannian metric on the parameter manifold M 
and thus associates a positive definite matrix j(t)g with each point 
in 9 £ M. The energy integral along the curve 7 is 



- f f Tr( 7 - 1 77" 1 7)) dF(6)dt. (4) 

' JO J M 

In this expression, dF(-) is a probability density on M, 7 means 
7(t)e, and 7 is the derivative of 7 with respect to t. Geodesies in Mi 
minimize _E 7 , and a variational approach is used in (3j to obtain the 



differential equation 7 

lit): 



77 7 for 7 (i), which implies 

7 (0)exp(7(0)- 1 7 (0)t) 



The right-hand side of this differential equation is observed to be a 
Christoffel symbol. 

The induced metric at a point a G S is 

G ff («,»)= / Tr( ff (a)-VWs(^rVM) 

where w and v are in the tangent space TS a of 5 at it and g* is the 
push-forward of g : S — > M. For a smooth curve 7 : [0, 1] — > S, 
the energy integral restricts to 



Enj 



Tr (ff(7W) _1 ff»(7W)ff(7W)" 1 3*(7(i))) dF(0)dt 



The geodesies, which are the extremal curves of E 1 , satisfy 

1 = - r 7(7,7) 

where r denotes the Christoffel symbol for the Levi-Civita connec- 
tion on S, 

In terms of local coordinates in S, geodesic equation in S may be 
obtained by solving a variational problem on the path u. To set this 
up, it is convenient to abuse notation and define a smooth function 
u : [0, l] 2 — > S such that u(s,t)\ s =o = u(t). With this notation, 
and coordinatizing G a as Qij, 



d_ 

ds 
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y~] Qi,j(u)u t u{ dt 



d k Qi,j(u)Usulul 



+ 2^Qi,j{u)u ts u 3 t dt 



If this expression is set to zero, further algebraic simplification leads 
to a differential equation (in coordinates) that characterizes geodesies 
in S: 



i,j \ k k 



* 3 

u t u J t . 



Returning to the example pictured in Fig. 1, the local coordinates 
in S — K 4 are x\, yi, X2, and 3/2- The positive definite matrix g(u) 
corresponds to the Fisher information matrix To given in (2). The 
inverses and derivatives needed are calculable, with 



R1R2 
kA(k) sin 2 (iy9i — if 2) 
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— sin ifij 
sin ifij cos (fij 



sin tfij cos (fij 

— COS 2 Ifij 



and 



where 



t-i \ ^ d . \ - d ^_ . 



d j- _ k,A{k) / sinifij 
dxj e ~ R A 3 \sin<pj -2cos^j 

d _ kA(k) 1— 2smifij cos ifij 
dy~ ~ ~RT V cos ^ 



Fig. [4] shows trajectories obtained for a particular case of the 
example scenario. The target is stationary at (1,1), and the sensors' 
prior distribution on the target location is normal with mean (1,1) and 
covariance 0.011. Sensor 1 starts at (0,1) and Sensor 2 starts at (1,0), 
and initial directions of motion are defined by the geodesic for this 
configuration. The sensors move in this direction for a fixed period 
of time, a new set of directions is determined from geodesic calcu- 
lations based on the new configuration, the sensors move again, and 
so forth. The dotted trajectories are extrapolations; they indicate the 
directions defined by the geodesic computation at the last iteration 
computed. 
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Fig. 3. Sensor trajectories based on geodesic approximation for the 
example scenario. Sensor 1 starts at (0,1) and Sensor 2 starts at (1,0). 
The target is at (1,1). 



5. DIVERGENCES ON M 

The proposed scheme for sensor management involves following, at 
least locally, geodesic curves in S defined by the Riemannian metric 
5* inherits from Ml. Geodesies maximize energy integrals of the 
form Q, so it is desirable to understand how optimization in this 
sense relates to the amount of information gathered by the sensor. 
Consider first the Kullback-Leibler divergence D(M2\ |A/i) for two 
multivariate normal distributions with equal means and respective 
non-singular covariance matrices g and h. This is given by 



lo« 



\h\ 
If I 



where | ■ | denotes determinant and I is the identity matrix. A diver- 
gence on A4 may be defined by 



Ao.(0,/») 



i TV (.,„- -,) + ilo 8 (|i| 



dF{0). 



Here, the two positive definite matrices g and h are regarded as aris- 
ing at each point of M from two Riemannian metrics. It is evident 
that zV (5, 3) = 0^nd d g A KL (g,h)\ h=g = d h A KL (g,h)\ h=g = 0. 
The corresponding Riemannian metric on A4 is 



'J 3 ^KL\ h=g 



d h An\h= 



Mg-gg-g) dF(0), 



as appears in (4]l. 

Similarly, one can define a divergence on M motivated by mu- 
tual information by 



A m {g,h) = I -{log 

' M 



20 + ff h) 



+ log( -(I + h^g) 



dF{9) 



This "symmetrized" mutual information expression is equivalent to 
1 



Am(g,h) = / \ log 

' M 



-(l + g~ l h) 



dF(6) 



As with A KL , it is clear that A m (g, g) — 0. Calculation reveals that 
d g A m (g, h)\h= a = d h A m (g, h)\h= g = and that the correspond- 
ing Riemannian metric on A4 is 



Thus, despite arising from different concepts of information (i.e., 
Akl from Fisher and Am from Shannon), both of these divergences 
give rise to exactly the Riemannian metric on Ml used in the geodesic 
computations of Section]?] 

6. CONCLUSION 

In this short paper, we have built upon results in differential geom- 
etry, outside the context of information geometry, to introduce an 
information-geometric approach to sensor management. The ap- 
proach begins with the observation that, when the goal of sensing 
is parameter estimation, the effect of selecting a particular sensor 
configuration amounts to imparting a Riemannian metric on the pa- 
rameter manifold M via the Fisher information. The collection of 
all such metrics is the Riemannian manifold M(M), for which the 
metric, geodesic equations, and other differential geometric aspects 
are known. With the assumption that our choices of sensor config- 
uration are parametrized by a smooth "sensor manifold" S, we ob- 
served that S inherits a Riemannian structure from M and used this 
to obtain a differential equation characterizing geodesic curves in S. 
In the purely geometrical work on which we have built, the measure 
on M is a volume form that corresponds to the statistical notion of a 
(minimally informative) Jeffreys prior. We observe that this may be 
replaced by an informative prior, as would typically be desirable in 



sensor management applications. Navigation along geodesic curves 
in a Riemannian manifold maximizes an energy integral involving 
the metric. We have constructed two distinct divergences on A4 cor- 
responding to familiar information-theoretic quantities (Kullback- 
Leibler divergence and mutual information) that have been used by 
various authors as criteria in designing sensor scheduling algorithms. 
Both of these are shown to lead to the same Riemannian metric on 
M, suggesting the information gathering merit of sensor schedul- 
ing based on following geodesic curves defined with respect to this 
metric. 

While the work presented here is mostly conceptual, we have 
shown enough specifics of how the proposed method manifests in a 
concrete example to indicate its feasibility. We are continuing to de- 
velop complete application examples while simultaneously working 
out rigorous specifics of some of the mathematical foundations. 
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